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Abstract: Inclusive momentum spectra and multiplicity distributions of charged particles measured with the 
BESII detector at center of mass energies of 2.2, 2.6, 3.0, 3.2, 4.6 and 4.8 GeV are presented. Values of the second 
binomial moment, R2, obtained from the multiplicity distributions are reported. These results are compared with 
both experimental data from high energy e + e~ , ep and pp experiments and QCD calculations. 
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1. Introduction 

Inclusive hadron production from e + e~ anni- 
hilations is very valuable for testing Quantum 
Chromodynamics (QCD). Perturbative QCD 
(pQCD) can give quantitative analytical predic- 
tions based on the Modified Leading Logarithmic 
Approximation (MLLA) £Q under the assump- 
tion of Local Parton Hadron Duality (LPHD) 0. 
Many experimental results have been reported 
from high energy experiments [3| @|, which are 
in good agreement with pQCD predictions. How- 
ever, hadron production has not been studied well 
at low energies due to insufficient data, particu- 
larly at the boundary region between pQCD and 
non-pQCD. It is therefore very interesting to test 
pQCD in this region with low energy e + e~ data. 

Here we study the inclusive momentum spec- 
trum, the hadron multiplicity distribution, and 
the second binomial moment. The inclusive mo- 
mentum spectrum is defined in terms of the vari- 
able £ = — ln(2p/y / s), where p and y/s are the 
momentum of the charged particles and center- 
of-mass (cm.) energy respectively. The sec- 
ond binomial moment a measure of the 
strength of hadron-hadron correlations and a 
sensitive probe for higher order QCD or non- 
perturbative effects [H], is defined as i?2 = 
(n c h(n c h ~ I)) /(ndi) , where n ch is the charged 
particle multiplicity. There has been a long 
standing discrepancy between the values of 
calculated to next to leading order (NLO) and 
those measured in e + e~ experiments [3]. 

In this paper, we report measurements of the 
inclusive momentum spectra and multiplicity dis- 
tributions of charged particles obtained with the 
upgraded Beijing Spectrometer (BESII) at the 
Beijing Electron-Positron Collider (BEPC) with 
center-of-mass energies of 2.2, 2.6, 3.0, 3.2, 4.6 
and 4.8 GeV. Results are compared with pQCD 
calculations, as well as those from high energy 
e + e~, ep and pp experiments. 

2. Detector and trigger 

The measurements were done using the data 
collected for the R scan |7I8| with BESII, a con- 
ventional solenoidal magnet detector that is de- 



scribed in detail in Ref. A vertex chamber 
comprising 12 tracking layers surrounds a beryl- 
lium beam pipe and provides input to the trigger 
system, as well as coordinate information. The 
primary tracking device is the cylindrical main 
drift chamber (MDC). It has 40 layers of sense 
wires and yields precise measurements of charged 
particle trajectories; it also provides dE/dx infor- 
mation which is used for charged particle identifi- 
cation. Outside the MDC, there is a barrel time- 
of-flight system (BTOF) consisting of an array of 
48 plastic scintillator counters that are read out 
at both ends by fine-mesh photomultiplier tubes. 
Electron and photon showers are detected in a 
sampling barrel shower counter (BSC) that cov- 
ers 80% of the total solid angle. It consists of 
24 layers of self-quenching streamer tubes inter- 
spersed with lead; each layer has 560 tubes. The 
outermost component of BESII is a muon identi- 
fication system consisting of three double layers 
of proportional tubes interspersed in the iron flux 
return of the magnet. 

The triggers and the measurement of the trig- 
ger efficiencies for Bhabha, dimuon and hadronic 
events are the same as those described in Ref. [7] . 
The efficiencies are measured by comparing the 
responses to different trigger requirements in spe- 
cial runs taken at the peak of the J/ip reso- 
nance, and are determined to be 99.96%, 99.33% 
and 99.76% for Bhabha, dimuon, and hadronic 
events, respectively. Their uncertainties are 
about ±0.5%. 

3. Hadronic event selection and back- 
ground subtraction 

The hadronic event samples are almost the 
same as those used for the measurement of R val- 
ues at these energies. Sources of background are 
cosmic rays, pair produced leptons, two-photon 
processes, and single-beam associated processes. 
First, clear Bhabha events are rejected from the 
sample. Next, hadronic events are selected. Spe- 
cial attention is paid to two-prong events, where 
cosmic ray and lepton pair backgrounds are espe- 
cially severe, and additional requirements are im- 
posed to provide extra background rejection [7J. 

The hadron selection proceeds by first selecting 
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charged tracks and then selecting events. For the 
track-level selection, we require the track polar 
angle satisfy | cos 6*| < 0.80, the track must not 
be identified as an electron or a muon, and the 
distances of closest approach to the beam in the 
transverse plane and along the beam axis should 
be less than 2.0 and 18 cm, respectively. The 
following criteria are also used to define "good" 
charged tracks: 

(i) p < E beam + (5 x op) , where p and E beam are 
the momenta of the track and the beam en- 
ergy, respectively, and a p is the momentum 
resolution for charged tracks at p = Ebeam ; 

(ii) E < 0.6Ebeam, where E is the energy in the 
BSC that is associated with the track; 

(iii) 2 < t < t p + (5 x at) (in ns), where t and 
t p are the time-of-flight for the track and 
a nominal time-of-flight calculated for the 
track assuming a proton hypothesis, and at 
is the BTOF time resolution. 

(iv) pt > 0.08GeV, where p t is the transverse 
momentum. 

The event-level selection requires at least two 
charged tracks with at least one "good" track 
satisfying the above requirements, the total de- 
posited energy in the BSC should be greater than 
Q.2&Eb eam , and the tracks selected should not all 
point along either the +z or — z direction. These 
criteria help reject background. For two-prong 
events, cosmic-ray and lepton pair events are re- 
moved by requiring that tracks should not be 
back-to-back and that there should be at least 
two isolated neutral clusters with E > 100 McV 
and with the differences in azimuthal angle with 
charged tracks more than 15°. This last require- 
ment rejects radiative Bhabha events. The above 
selection procedures remove most backgrounds. 
Backgrounds from two-photon processes are neg- 
ligible after hadron selection. 

Residual background contributions from 
Bhabha processes and tau pair production are 
subtracted using a Monte Carlo analysis, e.g. for 
Bhabha events, the number of Bhabhas in the 



hadronic events is 

pjpass 

Nbb = j^gen £ ■ Cbb, (1) 

where N^ b ass is the number of Bhabha events sur- 
viving hadronic event selection, N? b en is the total 
number of generated Bhabha events, C is the in- 
tegrated luminosity at each energy point, and abb 
is the total Bhabha cross-section calculated by 
the Monte Carlo program. A detailed compari- 
son between data and Monte Carlo data shows 
that the maximum discrepancy in any momen- 
tum bin is 10%; hence the error due to the Monte 
Carlo simulation of Bhabha events and tau de- 
cays is less than 1% in any bin, far less than the 
5-10% statistical uncertainty. 

The beam-associated background is treated by 
fitting the distribution of event vertices along the 
beam direction with a Gaussian for real hadronic 
events and a second order polynomial for the 
background, as described in Refs. [HlCOj- Tabled 
summarizes the background contributions to the 
hadron samples from Bhabhas, tau decays, and 
beam-associated events. The typical background 
is a few percent below rr threshold and about 
15% above. 

We do not subtract the contributions from the 
V>(4040), ^(4160), or the ^(4415) resonances to 
the high energy points. Using the resonance pa- 
rameters from the Particle Data Group [TTJ and 
varying them within errors, the contribution to 
the 4.6 GeV point is expected to be below 2.2 %, 
while that at the 4.8 GeV point is less than 1 %. 

Table 1 



Background contribution to the hadron sample. 





Nhad 


Bcam-assoc. 


e + e 


t pair 


(GeV) 




(%) 


(%) 


(%) 


2.2 


1410 


3.82 


0.61 




2.6 


4968 


3.72 


0.48 




3.0 


2030 


3.01 


0.56 




3.2 


1828 


4.53 


0.35 




4.6 


1315 


6.86 


0.18 


6.98 


4.8 


1282 


8.72 


0.14 


6.34 
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4. Hadronic event generator and initial 
state radiative correction 

Hadronic events are simulated by the JET- 
SET7.4 [T21 and LUARLW Q3| Monte Carlo pro- 
grams. LUARLW is used for energies at 3 GeV 
or lower. JETSET with tuned parameters is used 
in the region from 3 to 4 GeV, and JETSET 
with default parameters is used for 4 GeV and 
higher energy points [S]. The detector response 
is based on EGS for electromagnetic interactions, 
while for hadronic interactions, parameterizations 
are used. Above 3.77 GeV, the production of 
D, D*, D s , and D* is included in the genera- 
tor according to the Eichten Model ^31- A Monte 
Carlo event generator has been developed to han- 
dle decays of the resonances in the "radiative re- 
turn" processes e + e~ — > jJ/tp or jip(2S) [T5] . 
Four distributions comparing experimental data 
and Monte Carlo data at 2.2 GeV are shown in 
Fig. The data and Monte Carlo distributions 
agree reasonably well. 

Initial state radiation effects are simulated by 
a Monte Carlo program which is based on the 
scheme described in Ref. ^Hj- The error on the 
effective radiative correction factor (1 + 8 bs) is 
estimated to be less than 3% by comparing dif- 
ferent formalisms ^ZJEi anc ^ usm S different cuts. 



5. £ spectrum 

The charged particle inclusive hadron momen- 
tum spectrum can be expressed as a function of 
£ = — ln(2p/y / s). A purely analytical approach 
giving quantitative predictions for £ is a QCD cal- 
culation using MLLA under the assumption of 
LPHD expressed as 



1 da 

&had d£ 



= 2K LPHD x f M LLA(t,Aeff,N c ,n f ) (2) 



where Klphd is an overall normalization fac- 
tor describing hadronization, fuLLA is a complex 
function of £ and the effective scale parameter 
A-eff, N c is the color factor, and n/ is the num- 
ber of active quarks. Eq. |5] is only valid in the 
range < £ < ln(0.5yi/A e// ). 

The £ spectrum at each center-of-mass energy 
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Figure 1. Comparison between data (dots with 
error bars) and the LUARLW Monte Carlo (his- 
togram) at 2.2 GeV: (a) Multiplicity; (b) Spheric- 
ity; (c) Rapidity; and (d) Transverse momentum. 



is obtained using the following: 

1 dfJ -n(C\ 1 N °bs{£,i) /ox 

IF ~ ^^JTr a7 — ' W 

O had Clt,i iV had 

where, Nhad is the total number of hadronic 
events observed, -/V f, s (£j) is the number of 
charged tracks in bin-z, Af; is the bin width cho- 
sen commensurate with the detector resolution 
to avoid significant migration between bins, and 
C(£) is a multi-source correction factor includ- 
ing effects from the detector, ISR, and hadronic 
event selection. This correction is determined 
by two Monte Carlo samples: sample I (hadron 
level) includes neither initial state radiation (ISR) 
nor detector simulation and sample II (detector 
level) includes both ISR and detector simulation . 
Events of sample II are reconstructed in the same 
way as for data and subjected to the same selec- 
tion criteria. Distributions at the hadron level are 
only from particles with a lifetime greater than 
3 x 10 -10 s, which are considered as stable |19|. 

The correction factors C(£) for each bin are 
given by 



°VSiJ — I ATtotal )l\ iKTtotal I 



AJtotal 
gen 



N, 



(4) 



del 



where N*"^ 1 is the total number of generated 
Monte Carlo events for sample I, N d °* t al is the 
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number of events which pass the hadron selec- 
tion in II, and N gen (£i) and iVdet(£i) are the num- 
ber of entries in the ith bin, respectively. These 
correction factors are determined with the LU- 
ARLW and JETSET Monte Carlo simulations at 
3 GeV, where both Monte Carlo generators are 
valid. Their differences are included as system- 
atic errors for all center-of-mass energies. Other 
systematic errors include variations due to selec- 
tion criteria changes, etc. 

The measured £ spectra at the six energies be- 
tween 2.2 and 4.8 GeV are shown in Fig. For 
each point, we add the statistical and systematic 
errors in quadrature. These spectra are fitted 
with Eq. 0in the range 0.5 < £ < ln(0.5^/s/p ), 
where po ~ 0.35 GeV. The solid lines in the plots 
indicate the actual fitting range while the dashed 
line (outside the fitted range) indicate the QCD 
calculations using the fitted parameters, Klphd 
and A e f / , which are shown in Tabled At 2.2 GeV 
the fit range is very restricted and little of the 
peak region remains for comparison with MLLA; 
therefore this point is excluded. The fitted values 
of A e ff increase with decreasing E cm . The differ- 
ences from varying fitting ranges are included as 
systematic errors. 



Table 2 

The fitted A e // and Klphd values. The first 
uncertainty is statistical and second systematic. 



Ecm 


A eff (McV) 


Klphd 


2.6 


342 ± 7 ± 23 


1.523 ± 0.018 ± 0.023 


3.0 


325 ± 9 ± 25 


1.573 ± 0.027 ± 0.026 


3.2 


286 ± 1 7 ± 37 


1.532 ± 0.028 ± 0.052 


4.6 


239 ± 14 ± 32 


1.472 ± 0.029 ± 0.039 


4.8 


238 ± 1 5 ± 32 


1.482 ± 0.029 ± 0.038 
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Figure 2. Measured £ spectra (solid dots) at 2.2, 
2.6, 3.0, 3.2, 4.6 and 4.8 GeV. The solid curves 
show the data region fitted with the limiting spec- 
trum. The dotted line is an extrapolation of the 
fitted result. 



Eq. [21 has a maximum in each £ spectrum as 
shown in Fig. [21 Using the fitted parameters A e f f 
from Table |21 the peak positions, £*, can be de- 
termined directly from the limiting spectrum (Eq. 

EJ. 

The peak positions can also be obtained by 
fitting a Gaussian or distorted Gaussian [^Oj 
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over a limited range of 1/N event x dN/d^ > 
0.6(1/N event x dN/d£) max . We choose the Gaus- 
sian fit since the distorted Gaussian gives con- 
sistent results but with larger errors due to the 
larger number of free parameters. The changes 
found when the fitting range is varied are in- 
cluded as systematic uncertainties. Table 01 lists 
the fitted results. The results of the two methods 
(pQCD prediction of Eq. El and fitted Gaussian 
form) are consistent at higher energies, and the 
differences between the two methods are also in- 
cluded in the systematic uncertainties. 

Table 3 

Peak positions of the £ spectrum found by fit- 
ting with a Gaussian. The first uncertainties are 
statistical and the second systematic. 



Ecm (GeV) 


r 


2.2 


1.063 ± 0.017 ± 0.032 


2.6 


1.189 ± 0.012 ± 0.076 


3.0 


1.304 ± 0.017 ± 0.049 


3.2 


1.361 ± 0.021 ± 0.057 


4.6 


1.646 ± 0.025 ± 0.134 


4.8 


1.633 ± 0.037 ± 0.193 



The measured peak positions from this work 
are plotted in Fig. |3| together with those of high 
energy e + e~ and ep data. MLLA/LPHD predicts 
the energy dependence of the peak position, £*, 
as P] 

£* = 0.5Y + VcV- c, (5) 

where Y = ln(0.5s/s/A eff ) and c is 0.2915 
(0.3190) for three (four) active flavors. By fitting 
our data to Eq. , we obtain A e / / = 262 ± 9 
McV. This value is consistent with the results 
from OPAL EH, ZEUS E2 and CDF EHj, which 
are 263 ± 4, 251 ± 14(ep), and 256 ± I3(pp) MeV, 
respectively. However, our value disagrees with 
L3 El and the result of Ref. g|, which are 200 ±3 
and 232 ±3, respectively. Figure shows that 
is approximately linear in In ^/s. A straight line 
fit of the BES results to as a function of a/s 
gives a gradient of 0.779 ± 0.122. The gradients 



reported by OPAL [251, ZEU S ESI and HI E3 
are 0.637 ± 0.016, 0.650 ± 0.077, and 0.75 ± 0.05, 
respectively. Our result is somewhat larger. 



4.5 




Figure 3. Energy dependence of the peak position 
from Table [21 The curve is obtained fitting 
with only the BES data by Eq.0 but results from 
other e + e~ |3I lj and ep EH] experiments are also 
plotted. 



6. Momentum spectrum 

In order to study the behavior of low mo- 
mentum particles, MLLA/LPHD introduces a 
Lorentz invariant variable —^—E^-, where E 2 = 

o±ot dp 5 

P 2 + Ql and Q = 0.27 GeV EHj- The six BES 
Lorentz invariant, charged particle momentum 
spectra are plotted in Fig. QJ together with those 
measured by other experiments at higher energy, 
up to 130 GeV. The tendency of the spectra to ap- 
proximately converge at lower momenta is cited 
as evidence that hadron production at very small 
momentum p < 0.1 GeV is approximately energy 
independent. EH1- This behavior has been ex- 
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plained in Ref. |2E] as the coherent emission of 
low energy (i.e. long wavelength) gluons by the 
total color current. 
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Figure 4. Charged particle momentum spectra for 
center of mass energies from 2.2 to 130 GcV. Also 
shown are MARKI TASSO [3T], TOPAZ 
[32, and OPAL [55125] results. 



7. Charged particle multiplicity and the 
second binomial moment 

According to a NLO QCD calculation, the sec- 
ond binomial moment R2 is given by [5] 



R2 



11 



,-(v^) 



(6) 



with c = 0.55(0.56) for five (three) active flavors. 

The charged multiplicity distributions are ob- 
tained using a method which was first employed 
by TASSO |52| and modified by HI [27]. Three 
steps are used to correct the measured unnormal- 
ized charged multiplicity distribution Nfj!F(i), 



where i is the number of observed tracks. The 
first step is to correct for detector effects and 
selection criteria. It is performed using Monte 
Carlo sample II with its hadron level information. 
For each event, the number of observed tracks 
j with distribution N^f(j) is compared to the 
number of generated k tracks with distribution 
Ngfc (k). Denoting Njk as the number of events 
generated with k tracks when j tracks are ob- 
served and Nj as the number of observed events 
with j tracks, the correction matrix M(k,j) is 
expressed as 

M(k,j) = N jk /N j . (7) 
This matrix relates N^f(j) to N^(k) by 



•gen 



E 



M(k,j) 



N MC 



(J) 



(8) 



The second step is to correct for the pres- 
ence of QED initial state radiation which results 
in a reduction of the nominal cm. energy and 
thus changes the charged multiplicity. Therefore 
another set of correction factors is calculated ac- 
cording to 



Cf(&) = PNR(k)/p gen (k), 



(9) 



where PNn(k) is the normalized multiplicity dis- 
tribution for events generated at fixed cm. en- 
ergy without detector simulation and p gen (k) 
is the distribution N^^(k) after normalization. 
Finally, the corrected multiplicity distribution 
reads 



N™ r ; {k) = c F (k)Y,M{k, 



(10) 



The final step is done to reduce bias that 
might be introduced by the Monte Carlo gener- 
ator used [27J. An iterative procedure is used 
where the predicted multiplicity at the hadron 
level is reweighted using the previously unfolded 
multiplicity. This is repeated until convergence is 
reached. 

The measured multiplicity distributions at dif- 
ferent energies are shown in Fig. [5J along with the 
Monte Carlo predictions. Figure shows the en- 
ergy dependence of the mean multiplicity defined 
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as (n ch ) = n chP(n c h), where P{n ch ) is the 

n ch =0 

probability to have n c h tracks. 
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Figure 5. Multiplicity distributions at energies 
of 2.2, 2.6, 3.0, 3.2, 4.6 and 4.8 GeV. Solid dots 
are from this work; histograms are predicted by 
Monte Carlo. 



Based on the measured multiplicity distribu- 
tions, we obtain the second binomial moments 
i?2 , which are displayed in Fig. [7| together with 
both NLO calculations and published data at 
higher energies up to 100 GeV from other e + e~ 
experiments. Our measured R2 values, though 
with large errors, are consistent with other mea- 
surements at higher energies. The R2 value 
predicted by leading order QCD is significantly 
higher than the measured data, and while the 
NLO calculation comes closer to the data, the re- 
maining disagreement of about ~ 0.07 in R2 may 
indicate that R2 is a sensitive probe for higher 
order QCD or non-perturbative effects. 
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Figure 6. Energy dependence of mean multiplic- 
ity for energies below 8 GeV. Also shown are re- 
sults of ADONE [SI and MARKI [HOj. 
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Figure 7. Energy dependence of the second bi- 
nomial moment i?2- For other e + e~ data, see 
References |3I4| . 
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8. Summary 

Using BESII R scan data, we have measured 
the inclusive momentum spectra, multiplicity dis- 
tributions, and the second binomial moments at 
y/s ~ 2 - 5 GeV in the continuum. These results 
are compared with data from other experiments 
at higher energies, and with QCD model calcula- 
tions. 

The parameters Klphd and A e f f obtained by 
fitting the limiting spectra are compatible with 
those of the ZEUS, CDF, and OPAL experiments. 
However, the deviations become larger at energies 
below 3 GeV. 

The second binomial moment i?2 determined 
by the BES experiment, though with large errors, 
are lower than what is predicted to NLO. This is 
consistent with the long standing discrepancy be- 
tween the NLO calculation and high energy data 
from e + e~ experiments 3 . 

To further reduce the uncertainties of our mea- 
surements and test QCD model calculations to a 
precision of a few percent, higher statistics data 
collected with a better detector are needed. This 
could be one of the interesting physics goals of 
the CLEOc and BESIII projects |36I37| . 
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